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Abstract: We show that long-time, long-distance fluctuations of plane-symmetric hori- 
zons exhibit universal hydrodynamic behavior. By considering classical fluctuations around 
black-brane backgrounds, we find both diffusive and shear modes. The diffusion constant 
and the shear viscosity are given by simple formulas, in terms of metric components. For 
a given metric, the answers can be interpreted as corresponding kinetic coefficients in the 
holographically dual theory. For the near-extremal Dp, M2 and M5 branes, the computed 
kinetic coefficients coincide with the results of independent AdS/CFT calculations. In all 
the examples, the ratio of shear viscosity to entropy density is equal to h/ {AnkB), suggesting 
a special meaning of this value. 



Keywords: 



Black Holes, p-branes. Thermal Field Theory, AdS-CFT Correspondence. 



Contents 



|l]. Introduction ^ 

|2|. Diffusive and shear modes from horizon fluctuations ^ 
|2.1| Backgrounds ^ 



23 Current | 
|2.3| Boundary conditions ^ 
Quick derivation of Fick's law ^ 



|2.5| Checking the assumptions g 

|3|. Shear viscosity § 

Applications 11 

[4.1| Near-extremal D3-branes 11 



4.2 M2 branes 12 



4.3 M5 branes 13 



4.4 Dp branes 13 



4.5 A universal lower bound on r//s? 14 



H. Transport coefficients from AdS/CFT [16 
|5.1| R-charge diffusion 16 



5.2| Shear viscosity |18 

|6|. Discussion 19 

\A\ . Dimensional reduction 20 

0. Hydrodynamics 21 

in. Shear mode damping constant for the Buchel-Liu metric 22 



1. Introduction 

The holographic principle [|^, ^, ^ is believed to be fundamental to the construction of the 
quantum theory of gravity. It asserts that the physics in a region of space, in a theory 
with gravity, can be described by a set of degrees of freedom (some "dual theory", which 
does not contain gravity), associated with a hypersurface of smaller dimension in that 
space (the "holographic screen" ) . The nature of the degrees of freedom on the holographic 
screen is not specified by the holographic principle. In the explicit example of holography 
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discovered in string theory — the AdS / CFT duahty @, |5|, |6|, |^] , the holographic degrees of 
freedom form a local quantum field theory (A/" = 4 supersymmetric Yang-Mills theory in 
four dimensions), but in general this need not be the case. 

In this paper we elucidate the infrared properties of theories whose gravity duals con- 
tain a black brane with a nonzero Hawking temperature. We argue that the infrared 
behavior of these theories is governed by nothing other than hydrodynamics. Namely, by 
considering fluctuations around generic black brane solutions we found modes with dis- 
persion relations of the type u) = —iDq^, which are naturally interpreted as diffusion of 
conserved charges and damped hydrodynamic shear flow in the holographically dual the- 
ory. -"^ 

This result can be anticipated from the holographic principle. Indeed, the dual theory is 
at a finite temperature, equal to the Hawking temperature associated with the gravitational 
background, and is translationally invariant. On general grounds, we expect the infrared 
behavior of such a theory to be governed by hydrodynamics. The existence of the diffusive 
modes is consistent with this expectation. 

Some hydro dynamic-like properties of event horizons have been known for a long time 
in the framework of the "membrane paradigm" ^ . In this approach the stretched horizon 
is interpreted as a fluid with certain dissipative properties such as electrical conductivity 
and shear and bulk viscosities. Normally, the "membrane paradigm" is applied to black 
holes. However, for black holes the relation to hydrodynamics does not go beyond the level 
of analogy, in particular because their horizons lack translational invariance. The connec- 
tion to hydrodynamics is much more direct in the case of black branes. The membrane 
paradigm suggests the mapping between the bulk fields and the currents, whereas the long 
wavelength limit has to be taken to derive the hydrodynamic equations. 

Using the AdS/CFT correspondence, it has been shown |ll|, Q that the 

theories living on non-extremal D3, M2 and M5 branes behave hydrodynamically in the 
infrared. In this paper, we will see hydrodynamic behavior emerging as a common feature 
of black brane backgrounds. Our calculations do not rely on a particular realization of 
holography such as AdS/CFT. We find explicit formulas expressing the transport coeffi- 
cients (diffusion rate, shear viscosity) in terms of the components of the metric for a wide 
class of metrics. 

The paper is organized as follows. In section |2| we consider small fluctuations of 
black brane backgrounds, and demonstrate the emergence of Tick's law and diffusion from 
Maxwell's equations. We derive an explicit formula for the diffusion constant. In section ^ 
we extend the discussion to gravitational perturbations and compute the shear viscosity. 
In section Q the derived formulas are applied to various supergravity backgrounds. In the 
D3, M2, and M5 examples, the answers coincide with the kinetic coefficients in the dual 
theories, computed previously by AdS/CFT methods. In section ^ the kinetic coefficients 
in the theories dual to supergravity on Dp brane backgrounds are computed using the 
AdS/CFT recipe. The results coincide with those obtained in section 0. Section ^ discusses 
some unsolved problems and possible directions for future work. 

^The gravity dual of the sound wave, with the dispersion relation lo — VsQ — ijq'^, should also exist, but 
is beyond the scope of this paper. 
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2. Diffusive and shear modes from horizon fluctuations 
2.1 Backgrounds 

We start by considering linearized perturbations of a black-brane background. We take 
the Z)-dimensional background metric Gmn to be of the form 

p 

ds^ = Goo(r) dt" + Grr{r) dr^ + C(r) ^(dx')^ + Z{r) Kmn{y) dy"'dy'' , (2.1) 

1=1 

where the components Goo{r), Grrir), Gxxir), and the "warping factor" Z{r) depend only 
on the radial coordinate r.^ We assume that the metric has a plane-symmetric horizon 
at r = ro, that extends in p infinite spatial dimensions parametrized by the coordinates 
X*. In the vicinity of the horizon Gqo vanishes, Grr diverges, and Gxx, Z stay finite. The 
coordinates y"^ parametrize some d-dimensional compact space. 

We will be interested in small fluctuations of the black brane background, Gmn+Gmn, 
where Gmn is a small perturbation, and indices M, N run from to (1 + p + d). As 
usual, such perturbations can be divided up into "scalar", "vector", and "tensor" parts: 
Gmn, G^m, and G^u, correspondingly. Here indices fi, v run from to p + 1 (labeling the 
"non-compact" coordinates t, x* and r), and indices m, n run from 1 to d (labeling the 
"compact" coordinates y^). For our purposes we will be interested in vector and tensor 
perturbations only; this is motivated by the fact that effective hydrodynamics in field 
theories is constructed in terms of conserved currents and energy-momentum tensor of the 
theory. 

As in the Kaluza-Klein mechanism, the problem of analyzing y-independent ten- 
sor/vector perturbations reduces to the problem of gravitational/vector fields propagating 
on a lower-dimensional background. We write the metric of the dimensionally reduced 
background as 

p 

ds^ = gooir) dt^ + grr{r) dr^ + fc(r) Y^^dx'f . (2.2) 

Its components g^u are just G^u of ( |2.1D , multiplied by Z(r)'^/^. (Details of dimensional 
reduction are sketched in Appendix ^.) 

We will restrict ourselves to the linearized equations of motion for G^y and on 
the background ( |2.2D . The dynamics of the vector perturbations is governed by Maxwell's 
action 

W~ / d^+^x^ (^l.F^^-F^,^ , (2.3) 

where F^,^ = d^Ay — d^A^, with the gauge field A^^ being proportional to G^n- In (^), 
geflf is an effective coupling, whose overall normalization will not be important. If Z{r) 

^ Black brane metrics are solutions to the low-energy string- or M-theory (ten or eleven- dimensional 
supergravity) equations of motion. Full solutions also include other tensor fields such as the Ramond- 
Ramond field and/or the dilaton. However, to the linear order in gravitational perturbations, the presence 
of these extra fields does not affect subsequent discussion. 

^Our choice of metric excludes rotating branes. In other words, we do not consider nonzero chemical 
potentials. 
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is not constant (the compact space does not factorize), then g'^g is r-dependent, and is 
proportional to Z(r)~^P~^'^^^P . 



The dimensionally reduced metric (2^) has an event horizon at r = ro, in the vicinity 
of which goQ vanishes and Qrr diverges, 

goo = -{r - ro) 70 + 0{{r - r^f) , (2.4a) 

g,, = + 0(1) . (2.4b) 

r - ro 

Here 70 and 7^ are some positive constants. We assume that gxx{r) is a slowly varying 
function of r which does not vanish or diverge at the horizon, and is of the same order as 
goo when r — ro > ro. The Hawking temperature associated with the dimensionally reduced 
background (2.2) is equal to the Hawking temperature of the full metric (2.1), and is given 



by 

r = -^,/^. (2.5) 

47r Y 7,. 

We start by analyzing vector perturbations. Tensor perturbations Gfj,i, will be consid- 
ered in the next section. 

2.2 Current 

We shall illustrate the appearance of hydrodynamics on a simple problem of an Abelian 



gauge field on the black-brane background (2.2). Maxwell's equations, which follow from 
the action (|2.3D , read 

dJ^V^F'^'')=0. (2.6) 

To simplify the expressions, we will take g^s to be constant; thus Maxwell's equations 
become i\/~9 F^'^) = 0. The results for the r-dependent g'eff can be easily restored by 
replacing \/—g —>■ ^/—g/g1^ in the final answers. 

We will be considering fluctuations of the gauge field, which behave as oc e"*'^*"'"*'^' ' ^. 
Equations of motion for the field A^, together with appropriate boundary conditions give 
the dispersion law lo = uj{q). In the limit q — > 0, the dispersion law of the form co{q) = 
—iDq^, with D constant, is a sign of diffusion. We will be looking for the diffusive dispersion 



law on the black brane background (2.2). 

In principle, one can derive the dispersion law by directly solving Maxwell's equations. 
Here we will follow a somewhat different path. A conserved "current" can be defined 
directly in terms of the field strength One can show that if the gauge field satisfies 
the equations of motion with relevant boundary conditions, then at long distances Pick's 
law J* = —Ddij^ is valid. Then current conservation = implies that the charge 

density satisfies the diffusion equation dtj^ = Ddidij^. 

To define the current we need to introduce the stretched horizon |^. In our case, the 
stretched horizon is just a flat spacelike surface located at a constant r = rh slightly larger 
than rg: 

^'h > ''o , r^-ro <^ro. (2.7) 
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The outward normal to that surface is a spaceUke vector n^, which has only one nonzero 
component rir = gll'^iri^), so that n? = g^^nrUr = 1. The current which is associated with 
the stretched horizon is defined as in the original membrane paradigm approach ||8|, ^: 

r = n^F^r, ■ (2.8) 

The antisymmetry of F^^^ implies n^j^ = 0. That is, the current is parallel to the horizon. 
Further, Maxwell's equations, and the fact that all metric components g^y depend on r 
only, imply that is conserved, d^j^ = 0, for any choice of r^. The components of the 
current in the vicinity of the horizon are 

/ = F°V = Ynr--^,^ 

7o7r yrh-roYi'^ 

gxxlr 

The diffusion equation can be derived if we can prove Pick's law. We shall do that in two 
steps: 

• From the incoming-wave boundary conditions we infer a relation between electric and 
magnetic fields near the boundary: 

Fir = -.[^^^, r-ro^ro. (2.10) 
V 70 r - ro 

This relation is analogous to the relation B = — n x E for plane waves on a non- 
reflecting surface in classical electrodynamics. 

• We then show that if varies slowly in space and time, then F^i is proportional 
to diFor (analogously to E = — V(/9 in electrostatics), for appropriate choices of r^. 
Then the relations (2.9a), ( |2.9bD imply = —Ddij^, and give the expression for D. 



We now present details of the derivation. 
2.3 Boundary conditions 

Because of translation invariance of the metric ( |2.2| ), we can restrict ourselves to field 
configurations which are plane waves with respect to the spatial coordinates x*. Without 
loosing generality, we choose the wavevector to lie along the x = axis, so 

A^ = A^{t,r)e''^^. (2.11) 

In the hydrodynamics limit q is arbitrarily small. In particular we will assume it is much 
smaller than the Hawking temperature: g <C T. 

Fick's law j* = —D dij^ violates time reversal. Its origin thus must be rooted in 
the irreversible nature of the horizon. To incorporate time irreversibility we impose the 
incoming-wave boundary conditions on the horizon: waves can be absorbed by the horizon 
but cannot be emitted from there. 
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To derive the relation between electric and magnetic fields, we first write down the 
relevant Maxwell's equations and Bianchi identities, taking into account the form ( |2.11| ) 



V^gg'^^g^^dtFo. + d^V^ g'"^ g^^ Fr,) = , 
driV^g^'-g'^'For) + dA = , 

dtFrx + dxFor — drFox = . 



(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 



Using Eq. ( 2.12d ) one can write down field equations containing only the electric fields Fq^. 
and Fqx, 



dfFor - goo g^'^'q For - {iq)goo g^'^'drFo^ = , 
dtFox + goo dr(,g'~''drFox) - goo dr{g'''' d^F^r) = . 



(2.13a) 
(2.13b) 



In the near- horizon region r — tq ^ tq, these equations simplify considerably. Let us assume 
that fields vary over a typical time scale F"^, so that 3| ~ F^. At small (r — ro) Eq. ( p.l3a| ) 
implies 

(2.14) 



r-ro q 

J^Ot ~ ■^(Jrl'Ox ■ 



ro r2 

Hence at sufficiently small r — ro, namely, when 

r2 



ro 



(2.15) 



the third term in the left hand side of Eq. ( p.l3b| ) is small compared to the second term. 
This equation gives a "wave equation" for Fqx alone, 



dtFox - — (r - ro) dr[{r - ro)drFox\ = . 

7r 



The general solution to that equation is 



Fox{t,r) = h 



t + \ I — ln(r — ro) 

70 



+ /2 



ln(r — ro) 



(2.16) 



(2.17) 



where /i and /2 are arbitrary functions. The incoming-wave boundary condition picks up 
the /i term. This means 

Flo, 



dtFox = \ —{r - ro)drFox ■ 

V 7r 

From Eq. (|2.12d| ) and Eq. ( \2.l4) one finds that 

F)V Fox 



lor -rQ 



(2.18) 



(2.19) 



is independent of t. As we expect the solution to decay as t — > oo, this expression is zero. 
Thus we find Eq. (|1|). 
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2.4 Quick derivation of Fick's law 

From now on we will work in the radial gauge Ar = 0. We shall show that the parallel (to 
the horizon) electric field is dominated by the gradient of the scalar potential, Fqx ~ —dxAo, 
and that the normal (to the horizon) electric field is proportional to the scalar potential 
itself, Aq ~ For- Hence one finds Fqx ~ dxFor, and therefore ~ dxj^- 
When q = 0, Aq satisfies Poisson equation (see Eq. ( |2.12c| )) 



5r-(^/^/V°9,Ao) =0. 

With the boundary condition 74o(r) = at r ^ cxd, it can be solved to yield 

oo 

/ 9oo{r') grr{r') 



(2.20) 



Ao{r) = Co dr 



(2.21) 



Now let q be nonzero but small. Our first assumption, (to be checked later), is that unless 
r is exponentially close to vq one can expand in ^ series over q^ /T"^ <^ 1, and the leading 
term has the same r-dependence as in Eq. (2.21), with Cq depending on coordinates and 



Ao = A^^^+aS^ + 



4°)(i,x,r) = C7o(t)e^'?" I dr 



(2.22) 



When r ~ ro (but not exponentially close to rg) the ratio of FQr and is a constant which 
can be read off from the metric, 



For 



V-9{ro) 
gooiro) grr{ro) 



dr 



gooir)grrir) 



(2.23) 



Note that this ratio is finite since googrr is finite as r — > tq. Our second assumption is that 

\dtAx\ <C \dxAo\ . (2.24) 



Now Fick's law arises naturally: 





gxx^/lo{r - ro) 
dxAp 
gxx^/lo{r - ro) 

^OrJ gxx\/lo{r - ro) 



(2.25) 



-Ddxf , 
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where 



oo 



D = ^^^^ fdr Srrjr) ^2.26) 

ro 

For the gauge field with r-dependent coupUng gcs{f)-, the diffusion constant is 

oo 



gxx {tq ) gl^ {ro ) y^-gooiro) grr{ro) J \/-g{r) 



I -goo{r)grr{r)g,^{r) ^^.27) 
J \/-g(r) 



ro 



The natural estimate for the diffusion constant is D T^^. From Fick's law and continuity 
it follows that obeys the diffusion equation, and Co(t) oc e~'"*, F = Dq^. Notice that 
(2.15) requires that the stretched horizon has to be sufficiently close to the horizon, 



^«|^, (2.28) 



which is a stronger condition than rh — tq ^ tq. 



2.5 Checking the assumptions 



The short derivation of the diffusion law relies on two assumptions, Eqs. ( 2.22| ) and ( 2.24| ) 



We now verify that these assumptions are self-consistent. Namely, we assume Eq. ( 2.22| ) 
and then show that i) ( |2.24| ) is valid and ii) ( |2.22| ) satisfies field equations. 
From Eq. (|2.12a|) 

-rg^'A'o + iqg^^K = (2.29) 

(where prime denotes derivative with respect to r) one can find from Aq. Combining 
that with Eq. ( 2.22| ) and taking into account the boundary condition Ax\r=oo = 0, one 
finds 

a. = a(^^ + aW + ..., 

oo 

^(0) ^ Jlcoe-^t^^^^ [ dr' 9^^^'')9rr{r') (2.30) 

1 J ^FW) 



For r — ro ~ ro we have A^/Aq ~ F/(7. This comes from comparing Eqs. ( |2.30| ) and ( 2.22| ), 



taking into account that 500 ~ gxx for these values of r. However in the limit r — > ro the 
integral diverges logarithmically since g^r ~ — '^o)"^- For very small r — ro we then have 

^^~Ao-ln^^. (2.31) 

q r — ro 



Let us look at the condition ( p.24 ). We have 



' * —In — —. (2.32) 



\dxAo\ q^ r-ro 
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Since T^/g^ ~ /T"^ <C 1, the ratio is much smaller than 1 unless r — rg is exponentially 
small so that the logarithm is comparable to T^/g^. Therefore ( |2.24| ) holds if we choose 
the location of the stretched horizon rh so that 



In- 



,2 • 



(2.33) 



This means — cannot be too small. Still this condition can be satisfied simultaneously 



with (12:281) . 

Now let us check that Aq can be expanded as in Eq. ( 2.22 ). Expanding Eq. ( ^.12c ) in 
series over we find 



i(o) 



0. 



(2.34) 



Concentrating on r close to tq, and only on the orders of magnitude of the terms, we find 



q2aW ^ III. 



,^4°) + rqA(^' 



r -ro 

from which it follows that 



7rg 

r-ro 



q^ + In ■ 



A 



(0) 
' 



(2.35) 



.^In^ + r^ln^ "° 



r-ro 



r-ro 



A 



(0) 



ro r4 



A 





(0) 
' 



(2.36) 



where we have used the estimate gxx ~ 7o^o and F ~ q'^/T. We see that A^^^ is always 
smaller than A^^ for small q'^ /T'^. However, using Eq. ( |2.29[ ) one finds 



q yj ^ r — ro I r — ro ) 



(2.37) 



which is smaller than A^^x^ only if the condition ( ^.33 ) is satisfied. 

Thus, we have established that the assumptions underlying our derivation of Pick's 
law are valid, assuming that the location of the stretched horizon is chosen to satisfy ( 2.2^ ) 
and ( p3| ). 



3. Shear viscosity 

So far we have found that small fiuctuations of the stretched horizon have properties, which 
can be viewed as corresponding to diffusion of a conserved charge in simple fiuids. We now 
turn to the next simplest hydrodynamic mode - the shear mode. 

In principle, it should be possible to define the energy-momentum tensor T^"^ living 
on the stretched horizon in a manner similar to Eq. ( ^.g| ) and show that T^^ is conserved 
and, with a suitable choice for the location of the stretched horizon, satisfies the consti- 
tutive equations. We shall follow a simpler route. We will show that the corresponding 
fluctuations of the metric obey a diffusive dispersion law, to = —iDcf', and identify V with 
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rj/{€ + P) in the dual theory. Here rj is the shear viscosity, e and P are the equihbrium 
energy density and pressure (see Appendix We will call P the shear mode diffusion 
constant. 

We write ffuctuations of the {p + 2)-dimensional background ( p.2[ ) as g^u + and 
consider only those perturbations which depend on t, r, x = x^, with only two non- 
vanishing components of hf^i, 

hty = hty{t,x,r) , h^y = hxy{t,x,r) , (3.1) 

where x = x^, y = x'^ (assuming the black brane is extended along at least two spatial 
dimensions). This is motivated by the fact that hydrodynamic shear mode describes decay 
of the fluctuations in transverse momentum density of the fluid. We shall call perturbations 
of the type ( |3.1| ) gravitational shear perturbations. 

The linearized field equations for hty and hxy decouple from all other modes and so they 
can be treated separately."^ The easiest way to find the field equations for the gravitational 
shear perturbations is to notice that since hty and hxy do not depend on y, they can 
be viewed as zero modes of the Kaluza-Klein reduction on a circle along the y direction. 
Prom the standard formulas of Kaluza-Klein compactification one finds that the fields 
^0 = igxx)~^hty, Ax = {gxx)~^hxy satlsfy Maxwell's equations on a (p + l)-dimensional 
background, whose metric components are 

1 

9al3 = ga/siOxx)"-^ (3.2) 

The indices a, (3 run from to p. Indeed, the gravitational action contains the following 
piece 



^T^ia) ^ - 4 ^ Fc^pF^s i9xx)~' +.-., (3.3) 

where Fajj = daAp — dpAa, and the only non-zero components are Ftx-: Ftr, and Fxr- 

Thus the problem is reduced to that of Abelian vector field considered in the previous 

P 

section,^ with the identification l/ffg^ = {gxx)^^^- The shear mode damping constant V 
is given by the result ( p.27| ), evaluated with ^q^. Using the relation ( p.2p between g^^i, and 
(jf^i,, one finds 



V-9oo{ro)grr{ro) J gxx{r)^/-g{r) 



ro 



V-Goo{ro)Grr{ro) J Gxx (r) V^W) 

ro 

In writing the last expression, the relation between components of the black brane metric 
G^jy, and its dimensionally reduced version g^^i, was used, g^i, = G^yZ'^^^. 



*By a choice of gauge one can set hry to zero. 

^When solving for gravitational shear perturbations in 10 or 11-dimensional supergravity, one can con- 



sistently set fluctuations of all other flelds to zero. 
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Using thermodynamic relation e + P = Ts, one finds shear viscosity 

V-Goo(ro)G,,(i-o) J G„(r)y^GM 

This expression gives shear viscosity in terms of components of the black brane metric 
( |2.1| ), its Hawking temperature T and its entropy s per unit p-dimensional volume. 
The formula ( |3.5| ), and the corresponding one for the diffusion constant, 

oo 



^ _ y^GM Z(ro) r -Goo(r) G,,(r) 

G'.x(ro)x/-Goo(ro)G,r(ro) J ^/^G{r)Z{r) ' 

are the main results of this paper. 
4. Applications 

In this section, we shall apply the general formulas ( |3.5D and ( |3.6D to different gravitational 
backgrounds. For near-extremal D3, M2 and M5 branes, the results coincide with those 
found previously in Refs. [|lO|, 0, |l^ from the AdS/CFT correspondence. 

4.1 Near-extremal D3-branes 

The metric of a stack of N non-extremal D3 branes in type IIB supergravity is given in 
the near-horizon region by 

^2 o2 

ds^ = (-f{r) df + dx^ + dy"^ + dz^) + ^— - dr^ + R^d^l . (4.1) 

Here R is a constant, which depends on the number of D3 branes, R oc A^^/^, and 
/(r) = 1 — Tp/r^. (This metric is an important example, because for type IIB super- 
gravity on this background, the holographically dual theory is known explicitly.) The 
Hawking temperature of the background metric ( |4.1] ) is T = ro/vri?-^, and the entropy per 
unit (three-dimensional) volume is s = ^N'^T^. 

Applying Eqs. ( |3.6| ), ( |3.5| ), one finds the corresponding diffusion constant, 

T] = -N'^T^ . (4.3) 



and the shear viscosity 



The holographic dual theory for type IIB supergravity on the background ( |4.lD is 
= 4 supersymmetric Yang-Mills theory with gauge group SU{N), in the limit of large 
A^ and large 't Hooft coupling Q. This field theory lives in 3 + 1 (infinite, flat) dimen- 
sions, and must be taken at finite temperature, equal to the Hawking temperature of the 



gravitational background (4.1). Long-time, long-distance behavior in this theory is gov- 



erned by conventional hydrodynamics, with hydrodynamic variables being just densities 
®For systematic discussion of hydrodynamic fluctuations in supersymmetric theories, see | [l5| ]. 
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of conserved charges. This theory has 5'f7(4) global symmetry current (i?-current) which 
therefore must relax diffusively. The Pick's law for this current is ja = —Dab^j^, where 
a, h are adjoint SU{A) indices, which run from 1 to 15. When thermal equilibrium state re- 
spects the 5C/(4) symmetry (no chemical potentials for R-charges), the matrix of diffusion 
constants must be an invariant of the group: Dab = Dudab- Thus the R-charge diffusion 
is characterized by only one constant Dr, which can be computed by using the AdS/CFT 
approach [^]. The result of the AdS/CFT calculation of Dr is exactly equal to (^). 
Likewise, the result ( [4.3| ) exactly coincides with the shear viscosity of the supersymmetric 



Yang- Mills plasma, as computed by AdS/CFT methods llO|, 



In this example, the ratio of shear viscosity to entropy density is ry/s = l/(47r). 
4.2 M2 branes 

The metric of a stack of N non-extremal M2 branes in the eleven-dimensional supergravity 
is given by 

ds^ = H-^/^{-fdt^ + dx^) + H^/'^r^dr"^ + r^drj?) , (4.4) 

where H = 1 + R^/r^, f = 1 — r^/r^, and i? is a constant, which depends on the number 
of M2 branes, R oc N^^^. The Hawking temperature of the background metric ( |4.4| ) is 

Q 1 

(4.5) 



27rro H^'^iro) ' 
Using Eqs. (|3.6|) and (^), we find 



(4.6a) 



where 2^1(^1 b; c; z) is the hypergeometric function. In the near-extremal limit ro ^ R, the 
expression (l4.6a| ) gives^ D = 3/(47rr). 

The details of the CFT holographically dual to the near-horizon limit of the background 
( |4.4| ) are not known explicitly. However, it is known that its symmetry algebra includes 
(among other things) translation symmetry, and a global 5*0(8) i?-symmetry. Thus it 
makes sense to consider the long-time relaxation of charge densities of the corresponding 
conserved currents in the dual theory. The Minkowski AdS/CFT recipe [16|, |17] can be 



used to calculate correlation functions of the energy-momentum tensor and /^-currents in 
this theory at finite temperature, in the low- frequency limit. The AdS/CFT calculation |l2| 
shows that the i?-current relaxes diffusively, with the diffusion constant equal to 3/(47rT), 
and that the shear mode damping constant is l/47rT. Thus we see again that the general 



formulas (3^) and (|3.6|) reproduce previous AdS/CFT results. 



In this example, we also have tj/s = (47r) ^, independent of the ratio tq/R. 



^One can either use the fuU metric ( [4.4[ ) in the formula ( [j.6| ) and then take the near-extremal limit 
ro <^ R, or directly substitute the near- horizon limit of the metric {i.4) in Eq. (3.6), with the same result. 
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4.3 M5 branes 

The metric of a stack of N non-extremal M5 branes in the eleven-dimensional supergravity 
is given by 

ds^ = H-^/^{-fdt^ + dx^) + H'^/'\r^dr'^ + r'^d^l) , (4.7) 

where H = 1 + i?^/r^, / = 1 — r^/r^ , and -R is a constant, which depends on the number 
of M5 branes, R cx N^/^ . The Hawking temperature of the background metric (|4.7]) is 



47rro i?V2(ro) 
Applying Eqs. (3.6) and (|3.5|), we get 



(4. 



I^ = ^i^^/Vo)2F,(l,||-|), (4.9a) 

^ = ^^^/>o) = j^. (4.9b) 

In the near-extremal limit tq ^ i?, the expression (4.9a) gives D = 3/(87rr). Analogously 
to the previous example of M2 branes, these values reproduce the AdS/CFT results of 
Ref. |12]. Again, in this example rj/s = (47r)~-^. 

4.4 Dp branes 

Black p-brane metrics are solutions to the low-energy string theory equations of motion 
1 18]. The metric [p < 7) in the Einstein frame reads 

dsl = H-'-i^ [-fdf + dxl + --- + dxl) + {f-^dr^ + r^d^l_p) , (4.10) 

where 

r>7-p 7-p 

^ = 1 + ^' f = ^-3-p- (4-11) 



r 



The Ramond-Ramond field strength is given by 



R(7-p)/\ rl-^ + R7-P 
= {P- 7) hW^^p > (4.12) 



and the dilaton is 

The Hawking temperature is 



ij(3-p)/4(^) . (4.13) 



T=^—^H-'/\ro). (4.14) 
47rro 

According to the AdS/CFT dictionary, the isometrics of the (8 — p)-dimensional sphere in 
the metric correspond to a SO (9 — p) global R-symmetry of the dual theory. Fluctuations 
of the metric with one index along the brane and another one along the sphere become, 
upon dimensional reduction on S^~^, gauge fields propagating in the {p + 2)-dimensional 
bulk which couple to the R-symmetry currents on the {p + l)-dimensional flat boundary. 
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Dimensionally reducing the full Dp-brane metric ( [4.10| ) we obtain the following metric 

1 2(8-p) , o r, r.. 1,1 2(8-p) 

dsl^2 = Hvr~r- {-fdt^ + dxl + --- + dxl) + Hp^r^^ f'^dr^ . (4.15) 
The action of the SO (9 — p) gauge field is 

5~ [dP+^xy^^^2-^F^,, (4.16) 

where constant normalization factors are ignored. The effective gauge coupling constant 
depends on the radial coordinate, 

o P+1 _ 16 

5f^g(r) = H P r p . (4.17) 
Applying the formulas ( |2.27| ) and (^) we find 

We find again that rj/s = (47r)^^, and in the near-extremal regime vq <^ R the R-charge 
diffusion constant is D = {7 — p)/{8ttT). We shall reproduce these results in Section |5| by 
using the AdS/CFT prescription. 

4.5 A universal lower bound on 77/5? 

In all examples considered so far, there is a remarkable regularity: the ratio of shear viscos- 
ity to entropy density is always equal to (An)"^. This ratio is independent of the parameters 
of the metric and even of the dimensionality of space-time. This fact becomes more inter- 
esting if we notice that for any spacetime dimension the ratio r]/s has the dimension of the 
Planck constant. In the SI system of units the ratio r//s in all cases considered is 

- = — ^ « 6.08 X 10~^3 K-s. (4.19) 
s Anks 

One may suspect that the constancy of 7]/s is an inherent property of classical gravity, and 



hence Eq. (4.19) should be valid for any theory with gravitational dual description. This 
includes = 4 supersymmetric Yang-Mills theory in the regime of large A^, large 't Hooft 
coupling, and its deformations. 

We have checked this conclusion on the example of the supergravity solution recently 
constructed by Buchel and Liu |jl^ . It describes the finite-temperature J\f = 2* SU(A^) 
gauge theory at large 't Hooft coupling. The solution is a non-extremal generalization of 
the Pilch- Warner RG flow describing the mass deformation of the AA = 4 SYM by a 
hypermultiplet term. In addition to the temperature T, new parameters, the hypermul- 
tiplet masses m;,, mj, appear in the solution. The leading corrections to the metric are 
known in the high-temperature regime rrih/T <C 1, mf/T<^ 1, with rriij/T = 0, rrif/T = 
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corresponding to the near-extremal black three-brane solution. In Appendix |^ we show 
that the leading corrections to rj/s vanish identically. 

One may suspect that the number h/{4TTkB) is somehow special. We conjecture that 
it is a lower bound on r]/s. Since this lower bound does not contain the speed of light 
c, we suggest that this is a lower bound for all systems, including non-relativistic ones. 
This means we can check the conjecture on common substances where both viscosity and 
entropy density have been measured. 

Take, for example, water under normal conditions (298.15 K, atmospheric pressure). 



From the table [21| one finds s ~ 3.9 x 10^ J • K ^ • m ^ and t] ~ 0.89 mPa • s. This means 



rj/s 2.3 X 10^^'^ K • s, which is by a factor of 400 larger than ( 4.19 ). We have checked 
other common substances and found that the ratio i]/ s is always larger than h/{47rkB)- The 
minimum value we found is for liquid ^He at 1 MPa and 10 K, for which 7] ~ 3.07 (J.Pa • s 
and s ~ 4.95 x 10^ J • • m~^ [^], and the ratio rj/s is still larger than ( 4.1S| ) by a factor 
of 10. 

Let us now argue that in weakly coupled theories the ratio r]/s should be much larger 
than 1 and thus satisfies the lower bound. Weakly coupled systems can be described as 
dilute gases of weakly interacting quasiparticles. Finite-temperature A(/>^ theory with A ^ 1 
and gauge theories with small gauge coupling 5 <C 1 belong to this class. Nonrelativistic 
Bose and Fermi liquids at temperatures much below quantum degeneracy are also weakly 
coupled since they can be described as dilute gases of quasiparticles (phonons in the case 
of Bose liquids, and dressed fermionic quasiparticles in the case of Fermi liquids). 

The entropy density s of a weakly coupled system is proportional to the number density 
of quasiparticles n, 

s-^ksn. (4.20) 

The shear viscosity is proportional to the product of the energy density and the mean free 
time (time between collisions) r 

7?~ner, (4.21) 

where e is the average energy per particle (which is of the order of the temperature T). 
Therefore 



. (4.22) 

s Kb 

Now, in order for the quasiparticle picture to be valid, the width of the quasiparticles must 
be small compared to their energies, i.e., one should have 

- « e (4.23) 
r 

which means that 

^ » ^ . (4.24) 

s Hb 

Thus we see that (in units h = kB = 1) the ratio rj/s is large at weak coupling. Theories 
whose duals are described by supergravity are typically strongly coupled, thus naturally 
having r]/s of order 1. It is still puzzling that this ratio takes the same value of (47r)~^ for 
all such theories. 



-15- 



5. Transport coefficients from AdS/CFT 



We have seen that the two simple formulas (^^) and (|3.6| ) reproduce the known results for 
the transport coefficients of theories living on D3, M2 and M5 branes at finite temperature. 
For the general case of Dp brane there is no previous calculation to compare our result 
with. Here we show how to find the R-charge diffusion rate and the viscosity by directly 
calculating the correlation functions of the R-current and the stress-energy tensor following 



the prescription of Ref. [16|. (At zero temperature, such correlators were considered in 
Ref. E 



The metric of black Dp-branes is given by Eq. ( 4.10 ). In this section we work exclusively 
in the near-horizon limit r <^ R. It will be convenient to introduce a new radial variable 
(u = ro/r for metrics with even p and u = rg/r^ for the ones with odd p), in terms of 
which the near-horizon metric for even and odd p respectively becomes 

where f{u) = 1 — u^^^ (even p) and 

(7-p)(p+l) 

(5.2) 



(5.1) 



7-p 



where f{u) = 1 — u ^ (odd p). 

The dimensionally reduced metric ( [4.15 ), correspondingly, takes the form 



dsl = i?^ (^) {-fdt^ + dx? + • • • + d4) + R'-^^y^ {^y^"-' ^ (5.3) 
for even p and 

9—p 3 

dsl = i?^ (^) " i-fdt^ + dxl + .- + d4) + R^'^^y^ (^) " ^ (5.4) 



for odd p, with the same functions f{u) as in Eqs. (|5.lD , (|5.2|) . 

For the backgrounds (^]^), (|5.2| ) we now compute the R-current retarded correlators in 
dual gauge theories, as well as the correlators of the shear-mode components of the stress- 
energy tensor. For frequencies and momenta much smaller than the temperature we find 
that each of these correlators has a diffusion-type pole with a specific value of the diffusion 
constant. The calculations are very similar to those done for the D3-brane background in 
Ref. [^] and for M-branes in Ref. ||T^, where details of the approach can be found. Our 
results are summarized in Table ^. 

5.1 R-charge diffusion 

In computing two-point functions of R-currents in the large N limit, it is sufficient to treat 
the bulk gauge fields as Abelian ones, ignoring the self-interactions |24]. Accordingly, to 
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1 
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3 
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Table 1: Entropy density and diffusion constants for black p- and M-branes. There is no 
gravitational shear mode for Dl brane. For D5 brane, the entropy is linear in energy. 



find the R-charge diffusion mode, we consider Maxwell's equations 

1 
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(5.5) 



on the backgrounds (^!^), (|5.4D . Here to simplify notations we denote ^/—g = ^J^^g^yzj g\^. 
We choose the gauge = and restrict ourselves to fields which depend only on the u, t 



and X 



Afj,{u,t,z) 



Then for the components At and A^ one finds the system of equations 
g''u;A't-qg^^A', = 0, 

du ( - {u^qAt + oj^A,) = , 



(5.6) 



(5.7a) 
(5.7b) 
(5.7c) 



where prime denotes the derivative with respect to u. The equations for other components 
of Af^ decouple, and thus these components can be consistently set to zero. Eqs. ( ^.7aD and 
( 5.7b| ) can be combined to give a closed- form equation for A'^ 



d_ 

du 



-gg^^g^^ 



(5.8) 



According to the method developed in Ref. |11], one solves Eq. (5^) for A[ with the 
incoming- wave boundary condition at n = 1. This can be done in the hydrodynamic 
regime u/T <C 1, q/T <C 1. The obtained A'^ depends on a single overall normalization 
constant which can be fixed from Eq. ( |5.7b| ) and the boundary values A!^ = lim^^^o ^i(^)) 
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= lim^^o A.z{u). One then finds A'^ from Eq. ( 5.7a| ). Tliis is sufficient for computing the 
on-shell classical action, which is proportional to AtA[ — AxA'^ at li = 0. Taking derivatives 
of the action with respect to A^ and A^ one then obtains the current correlators. 
In backgrounds (^), ( ^.41) , Eq. ( ^ ) becomes 



d 

du 



(1 _ u° A 
du\ 



where for p even 

^even — 7 p , Oi^ygn — P 2 , Peven — 1 ; (5.10) 



while for p odd 



7 — p p — 1 

i^odd = , aodd = , Podd = . (5.11) 



The dimensionless parameters t», q are defined as 

Near the horizon, two local independent solutions of Eq. ( ^.9| ) are A[ ^ {u — 1)°"^ , where 
a± = ibito/i/. The "incoming wave" boundary condition corresponds to choosing a_ as the 
correct exponent. Solving Eq. ( ^.9| ) perturbatively in jb, q we find 

A't=C{u-l)-'^u->^ (l + wF'^P\u)+q'^G^'P^{u) + ---) , (5.13) 



here F^P\u) and G^P\u) are explicitly known (but rather cumbersome) functions, inde- 
pendent of to, q. The constant C can be found from Eq. ( ^.7b| ) and the boundary values 
of fields at u = 0. We have 

Ce^ «,+ q2 + •••)=(«, q ^0 + q2 A?) , (5.14) 

where the ellipses represent terms of higher order in to, q. Here we first see the emergence 
of the diffusion pole: C ~ {ivo + ^^q^)^^. This pole will appear in the correlators. It 
corresponds to the diffusion constant 



D = ^-^. (5.15) 



5.2 Shear viscosity 



The low-energy string theory equations of motion which give rise to the black p-brane 
backgrounds ( [4.10 ) can be deduced from the following Einstein- frame effective action 

S = j d'\V^ (^n - ^d^^dt^^ - ^e'^^i^n) • (5.16) 

Here ^ is the dilaton, Fn is the Ramond-Ramond field strength form (or its dual), and 
values of n and a are related to p for a given p-brane background. The equations of motion 
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that follow from Eq. ( 5.16| ) are 



= \ d,<^ d,^ + ^^^^ e'^* [f,„F- - g,^j , (5.17a) 

V^(e'^*F^-) =0, (5.17b) 



a 



□ $ = ^e''*F2. (5.17c) 
2n! 

To find the viscosity, we can substitute the perturbed metric gfj_u + h^p into Eq. ( ^.17a| ) , 
solve the resulting equations, and then use the Minkowski AdS/CFT correspondence to 
find the retarded correlators of the appropriate components of the stress-energy tensor in 



the dual gauge theory. The case of p = 3 was treated in this way in Ref . 1 11 ] . 

Alternatively, we may use the fact that, as noted in section the shear gravitational 
perturbations satisfy Maxwell equations with a coordinate-dependent effective coupling in 
a dimensionally reduced background. Thus the problem is reduced to the one treated in 
Section 5J: the system of the effective Maxwell's equations is given by (|5.7aD ~( ^.7c) ), with 



the metrics (5J), (|5.2|) for even and odd p respectively, ^/—g replaced by ^—gjg, 



eff 



gxxV~9 ~ (gxx) ^hty, Ax = (gxx) ^h-xy Thcu the equation for the effective A[ 

still has the form ( ^.9[ ), but with a new set of parameters, 

^even — 7 p , OL^^^n — 3 , Peven — P 6 , (5.18) 

'Y p p ^ 

i^odd = , Olodd = 2 , (3odd = ■ (5.19) 

The subsequent calculations follow the Maxwell case verbatim, and in the end we find 
the hydrodynamic pole in shear part of the correlator of T^^ at t» = —i<^jv -|- • • • . The 
dispersion relation is a; = —iVq^ + where V = I/AttT is independent of p, hence 
r}/s = 1/4tt. 



6. Discussion 

Let us briefly summarize. The analysis of small fluctuations around a black-brane back- 
ground shows that in the low frequency, long-wavelength limit, there exist perturbations 
that correspond to diffusive hydrodynamic processes on the stretched horizon. This re- 
sult indicates that the degrees of freedom holographically dual to a gravitational theory 
on a black-brane background (whatever their microscopic nature is) must have a rather 
conventional hydrodynamic limit at large distances and long times. 

We have also derived general expressions for the diffusion constant and shear viscosity 
in terms of the components of the background metric. For near-extremal Dp and M- 
branes, these formulas reproduce the results of the direct AdS/CFT calculations.* For all 

*It should be stressed that the notion of diffusion (more generally, hydrodynamics) on stretched horizons 
is not limited to the particular string-theoretic realization of holography. It is encoded in the low-frequency 
part of the spectrum of quasinormal modes for a given gravitational background (see for a discussion 
of quasinormal modes in this context). 
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non-extremal supergravity backgrounds considered in the paper we found that the ratio of 
shear viscosity to entropy density is a universal number^ equal to {4it)~^. 

We find these results rather intriguing, and some future work is yet to be done. First, 
at the moment we lack crisp understanding of why our analysis of small perturbations in 
the vicinity of the horizon yields the same result as the AdS/CFT calculation. 

Another question to be addressed is the complete treatment of gravitational pertur- 
bations. Indeed, the fluctuations of the stretched horizon must encode the full set of 
hydrodynamic modes, at least in the known AdS/CFT examples. This means that meth- 
ods similar to those used in Section |2| should reproduce propagating modes (sound), in 
addition to the non-propagating shear and diffusive modes. In particular, it is reasonable 
to expect that formulas analogous to ( p. 51 ) exist also for bulk viscosity and for the speed of 
sound. 

It would be interesting to extend and generalize our calculations to include other 
nontrivial supergravity backgrounds, such as rotating branes [26|, non-extremal Klebanov- 
Tseytlin solution [27|, or the full (as opposed to the high-temperature limit) Buchel-Liu 
metric |19]. 

Another interesting question is the gravitational description of non-liner terms in hy- 
drodynamic equations. Namely, the hydrodynamic constitutive relations for conserved 
currents contain in addition to the linear term (Tick's law), a non- linear convective term: 
J* = —Dd^j^ + (e + P)~^ j'^TT*. Here vr* is momentum density in the i-th direction; the 
last term is just charge density times velocity of the fluid. This term, which is quadratic 
in small fluctuations, is a direct consequence of Lorentz (or Galilean) invariance of the 
microscopic theory. The gravitational manifestation of the convective term may not be 
straightforward: it was shown in |15| that such non-linear terms give rise to 0{1/N'^) effects 
in correlation functions of conserved currents. In the context of AdS/CFT, this implies 
that gravity loop effects must be included. 

Finally, we hope that analyzing the hydrodynamic regime of black branes will lead to 
a better conceptual understanding of holography. 



A. Dimensional reduction 

In this Appendix we write down the formulas of the dimensional reduction scheme relevant 
for our discussion in the main text. More details can be found in |2|, |3^, 0. 



The dimensional reduction of a D-dimensional pure Einstein theory is performed by 
using the ansatz 

ds' = e'-Usj, + e'''^dsi., (A.l) 



^This observation traces back to the observation in [|l^ that shear viscosities and entropies of M-branes 
have the same A'' dependence. Curiously, the viscosity to entropy ratio is also equal to 1 /An in the pre- 
holographic "membrane paradigm" hydrodynamics there, for a four-dimensional Schwarzschild black 

hole one licLS ?7m.p. — l/167rGjv, while the Bekenstein-Hawking entropy is s = l/4Gjv. 
Sound waves from the AdS/CFT perspective were considered in 
Analogous non-linear terms also exist in the constitutive relations for stress T'-'. 
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where X is the lower-dimensional space and Y is the internal compactifying space with 
dimensions dx and dy respectively {dx + dy = D). In the reduced action, the lower- 
dimensional Einstein-Hilbert term and the kinetic term for the scalar (p will both have 
canonical normalization, if we take 



dy , I dx-2 ^^2) 



2{dx -2){dx + dy -2) ' V 2dy {dx + dy - 2) 



Considering vector- like fluctuations of the metric ( A. ID (gravitons with one index along 



the X space and one index along the Y one), one observes that the normalization of the 
corresponding gauge fields is not canonical. Schematically, 

I ^ I ./^e^(b-a)^ i^^l + • • • , (A.3) 

where constant prefactors are ignored. Thus in general the effective gauge coupling ~ 
^2(b-a)(j> ^^Yi be position-dependent. 

B. Hydrodynamics 

To make our presentation self-contained, here we briefly review the basic properties of 
hydrodynamic fluctuations. Hydrodynamics is an effective theory, which describes the 
dynamics of a thermal system on length and time scales which are large compared to any 
relevant microscopic scale. The degrees of freedom entering this theory, in the simplest 
cases, are the densities of conserved charges. 

As an example, consider a translationally invariant theory in flat space, which has a 
conserved current j^, which corresponds to some global symmetry. The hydrodynamic 
degrees of freedom are the charge densities j^, e = T^^, and vr* = T^^, where T'^^ is the 
energy-momentum tensor. The currents satisfy the conservation laws: 

dtf = -dif , (B.l) 
dte = -dy , (B.2) 
Oy = -djT^ . (B.3) 

The constitutive relations express the spatial currents j*, T*-' in terms of j'', e, and vr*. To 
linear order 

ji = -Ddd\ (B.4) 
Tij = 6ij {P + vide) - ^ ^ p 5ij dkiTk - [di-nj + dj-Ki - -5ij dkiTk) ■ (B.5) 

Here e = (e), P = ^{Ta) are equilibrium energy density and pressure, p is the number 
of spatial dimensions, Vg = {dP/deY/'^ is the speed of sound, and 5e = (e — e) is the 
fluctuation in energy density. The unknown kinetic coefficients D, ( and rj are diffusion 
constant, bulk and shear viscosities, respectively. The constitutive relation (|B.4| ) is often 
called Tick's law. 



- 21 - 



The linearized constitutive relations together with conservation laws give the disper- 
sion relations for hydrodynamic modes. Taking all charge densities to be proportional to 
^-lujt+iq-x^ one finds: i) diffusive mode, uj = —iq^D, ii) shear mode, uj = —iq'^r]/{e + P), 
and in) sound mode, u = v^q — iq^{C + rj)/{e + P). 

Finally, linear response theory allows one to extract kinetic coefficients from singulari- 
ties in equilibrium correlation functions of the corresponding conserved currents. This was 
used to find both diffusion constant and shear viscosity in some strongly-coupled theories 
from AdS/CFT recipe for correlation functions 0, 12]. We use this method in section 



C. Shear mode damping constant for the Buchel-Liu metric 

Non-extremal deformation of the Pilch- Warner RG flow recently found by Buchel and Liu 
1 19] is a gravity dual to the finite-temperature J\f = 2* SU{N) gauge theory at large N 
and large 't Hooft coupling. The five-dimensional metric is of the form 

where r is the radial coordinate in five dimensions. Functions A{r), B[r) satisfy supergrav- 
ity equations of motion governing the flow. The system of equations (which also involves 
two r-dependent scalars) is given explicitly by Eqs. (3.17) of ]|l^. Following ]jl9[ we intro- 
duce a new radial variable, y = e-^, y S [0, 1], with y = corresponding to the position of 
the horizon. The metric becomes 



ds^ = (-y2dt2 + dx") + dr^ ' . (C.2) 



Using Eqs. (3.17) of ]19] one finds the following expression involving the Jacobian dy/d' 



where prime denotes the derivative with respect to y. Functions p{y) and xiv) the two 
scalars satisfying supergravity equations of motion, and V is their potential. 



l/'l \'^1 1/^1 ^ 

V = — (-^ - p'' cosh2xj +— p«sinh2 2x-- f^ + yCosh2x) . (C.4) 



Applied to metric ( |C.2| ), our formula ( p.4[) for the shear mode damping constant reads 



Knowing A{y), p{y) and xiu)) one can in principle determine P using Eqs. ( p.4t) and 



( C.5 ). Unfortunately, the system of second order differential equations involving those 



functions is too complicated. Buchel and Liu ]19] were able to solve it perturbatively using 



ai (X (rrih/T)'^ <C 1 and 02 oc mj/T <C 1, where mi, and nif are masses of the bosonic and 
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fermionic components of the N = 2 hypermultiplet, as small parameters. To the leading 
order in ai, 02 the solution is 



A{y) = a - J ln(l - y^) + a?^i(y) + ai^2(y) , 
P{y) = 1 + ai/)i(y) , 
X{y) = a2X2{y) , 

where the scalars pi{y) and X2(y) obey linear differential equations 

(i-y^ (yx2)' + ^yx2 = 0, 

and functions Ai{y), A2{y) can be found by solving 



(C.6) 



y(l - y') A'l - (1 + 3y') A[ + 4y(l - y^) (p^ 



0, 



y(l - y2) ^'2' - (1 + 3y2) + ^y(l - y') (X2)' = . 

The scalars pi{y) and X2(y) are expressed in terms of the hypergeometric function 

,1 1 



Pi = {i-y')'/' 2F,{^,-,i;y') , 

X2 = {l-yyUF,{\,\,l-y^), 



(C.7) 
(C.8) 

(C.9) 
(C.IO) 

(C.ll) 
(C.12) 



and functions Ai, A2 are given by 



^1=6-4 



y zdz 



(1-^^)' 
4 P zdz 



(l-z2)2 



8-7r^ 
27r2 

8-37r 
Svr 



+ 



+ 



dx 



dx 



dpi \ (1 — X 
dx 



2\2 



dx 



(C.13) 



To the leading order in ai, 02, integration constants a, ^2 are related to the Hawking 
temperature of the black brane background by 



27r 



16 



1 + af Ui + ^ + "2 U2 + 



4 \ 



Svr / 



(C.14) 



For ai = 0, 02 = one recovers the original non-extremal black three-brane metric 

dy' 



ds^ = (2vrr)2(l - y2)-l/2 (_y2^^2 ^ ^^2) ^ 



2^2 ■ 



(C.15) 



Expanding the Jacobian to the leading order in ai, 02 we find 



dy 
dr 



ii-y') 
1 



l + ai(2pi + 2{l-yr{p[f 



^ 2\2/ / \2 1 y^ 



+ cxi[-xi + ^{i-yy{x2 



(C.I6) 
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Correspondingly, the shear mode damping constant ( C.5 ) can be written as a series expan- 



sion 



1 

V = ^ j dyy{l + alFiiy) + a^FaCy) + •••), (C.17) 



where ellipses stand for terms of higher order in ai, 02, and 



y 



4-4i2-^xi-|(l-2/m)^ + :f 
y 2 6 6tt 



(C.18) 
(C.19) 



where Ai = Ai{Ci = 0), = ^2(6 = 0). Using equations of motion (IOTD - (|Cj), (Q- 
( C.10| ) we now show that Fi and F2 vanish identically. 

To prove that Fi = 0, first we integrate Eq. ( C.9D from to y (taking into account 
that Ai{0) = 0, A[{0) = 0) to obtain the identity 



A^{y) = ^y{l-y')A[ + 2 J z{l - z'M] 



Inserting this into Eq. ( p.l^ ) we find 



(C.20) 



l-pi-{l-y')\p[r 



X-z" 



{p'lfdz. 



(C.21) 



To see that the right hand side of Eq. ( |C.21 ) is identically zero, we multiply both sides of 

by f 



the differential equation (C.7) by p'l- We have then 

yiplY 



(l-y2)2 2 



(C.22) 



This allows us to express the integrand in Eq. (C.21) as 

;{ply - {I - z')[{p[f]' 



1-^2 

-(ply-(l-z^f[(p[f]' + 4zil-zW,) 
2z^ 



l-z- 



^(pIY - 2zHl - z')[ip[f]' - izil - z'Mf 



Eq. ( C.22| ) ensures that the last line in Eq. ( |C.23D adds up to zero. Noting that 

-ipir - (1 - z^ufr + HI - ^')ip'if = ^ - (1 - zrip'if) 

and remembering that pi(0) = 1, p'i(O) = (see ( p. 11 )), we arrive at 

y 



1 -z^ 



ip[rdz = i-pi-{i-yr{p[f 



(C.23) 



(C.24) 



(C.25) 
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Consequently, Fi = 0. The proof of -F2 = is similar. We conclude therefore that (at 
least) to the next to leading order in the high temperature expansion parameters nib/T, 
m f /T the shear mode damping constant of the N = 2* gauge theory at large N and large 
't Hooft coupling is independent of the parameters of deformation (hypermultiplet masses) , 
and equals l/47rT. 
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